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"Method of searching for an object within a space, method and system of 
vectorial cartography incorporating this search technique, electronic apparatus 
employing this method of vectorial cartography, and a medium for vectorial 
cartography data obtained by the said method" 

5 

The present invention relates to a method of searching for an object within 
a space, as well as a search system implementing this method. It also relates to 
a method and a system of vectorial cartography incorporating this search 
technique, as well as electronic devices employing this method of vectorial 
1 0 cartography and media for vectorial cartography data obtained using the said 
method. 

In vectorial cartography, the objects present in a domain are described in a 
matrix which may be of a considerable size and whose processing may be 
particularly time-consuming. Moreover, when it is a matter of searching for 
15 objects recorded in a subdomain of the principal domain, it becomes necessary 
to construct an extract from the principal matrix. 

More generally, the problem posed is the determination of a subset of 
objects contained in a subdomain of a domain containing a set of semantically 
homogeneous objects of variable geometry, and the construction of a matrix 
20 describing the objects belonging to this subset. 

Consider a set Q of N objects that are semantically homogeneous and of 
variable geometry (shape and size). They are distributed randomly in the space 
OXYZ, where they occupy a domain A inscribed within a parallelepiped between 
(X min , Y min , Zmin) and (X max , Y max , Z max ), this parallelepiped constituting the limits of 
25 the "space" A occupied by the domain. 

In a two-dimensional system, the domain could be denoted by the term 
"territory", inscribed in a rectangle between (X min , Y min ) and (X max , Y max ). 



The N objects are described by a matrix M, which gives the shape for each one 
by means of a list of points (with coordinates) and topological descriptives 
(edges, faces, etc.). At the same time, the coordinates give the position of the 
object in space. Each object therefore has its own elementary space X 
represented by a circumscribed parallelepiped. The following observations can of 
course be made: 

- the elementary occupied spaces can intersect, 

- (Xmin, Y min , Z min ) and (X max , Y max , Z max ) are the minimum and the 
maximum respectively of at least one elementary parallelepiped. This is true if 
there is no geometrically explicit definition of domain A. 

It is required to determine the subset © of the v objects contained in a 
subdomain 6 inscribed in a parallelepiped between (x min , y m in, z min ) and (x maX) ymax, 
Zmax) inside the parallelepiped circumscribing A, and to construct the matrix \i 
describing them, the said matrix being an "extract" of M. 

According to the problem posed, an object contained in 5 can be regarded 
as meaning an object that is included completely or partially (object "touched" by 
or intersecting 8). 

Existing methods of speeding up the process will now be described: 

If the order of description of the objects of Q in the matrix M is random, to 
construct u. it is necessary to pass through all the N objects so as to verify for 
each one whether it belongs to 5 and if so, increment each time v. 

To speed up the construction of the extract ja, it is necessary to transform 

M. 

A first existing method of acceleration is preliminary sorting. Preliminary 
sorting of the objects of M makes it possible to speed up the search. The method 
used most widely is ordering of the objects in relation to the point (x min , y mi n, z min ) 
of the parallelepiped of their elementary space X. This makes it possible to 
reduce the number of complete checks of membership, by immediately 
eliminating all the objects for which this point is above (x max , ymax, 
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Zmax) of 8; the process can be applied a second time for the other end. 
Membership is then checked for the remaining objects. There are several 
variants of this approach. 

It will be noted that even if the complexity of certain checks is reduced (the 
5 smaller that 8 is relative to A, the greater is the reduction), it is necessary to go 
through all the N objects of M. Its avoidance would require the use of a 
dichotomizing search. 

A second type of method of speeding up the construction of the extract \i 
consists of breaking down the domain A into £ subdomains Aj called tiles, hence 
1 0 the name tiling. Objects straddling one or more tiles (subdomains) will be cut at 
the boundaries and additional information (regarding connection, to ensure 
semantic continuity) will be recorded. The original matrix M will therefore be 
broken down into S matrices Mj each having Ni objects. Obviously, 

15 2 Nj > N 

The domain 8 searched will cover ct tiles Aj. The efficiency of the method 
will depend on the number S of tiles of the original domain A and on the size ratio 
between 8 and the Aj, and therefore between a and S. On the other hand, 
20 unfortunately, the more we increase S, the more we increase the number of 

objects intersected and therefore the size of M. At the end, it will be possible for 
any objects that have been intersected to be reconstituted, thus making the tiling 
transparent. 

Tiling immediately seems more beneficial than simple sorting, especially 
25 as the latter can be applied to each tile. This is why it is widely used, for example 
in cartography. 

However, it has a major drawback in that it is necessary to cut the objects, 
and hence increase the total number of objects processed, consequently 
requiring the creation and use of boundary information, all of which complicates 
30 the processing algorithms. 



The proposed objective, relative to the state of the art, is to speed up the 
process for obtaining the extract without altering the objects in any way, so as not 
to increase their number and hence the size of the matrix M containing, for each 
object, a list of points and topological descriptives. 

To achieve this objective, the information on the elementary space 
occupied by each object is taken into account, this information either being 
contained initially in matrix M, or it must be calculated and then recorded in the 
matrix, and the set of objects is passed through a number of "sieves" of finer and 
finer mesh, so as to obtain groups of objects of particular sizes. 

Thus, a method is proposed for searching for an object contained in a 
domain 8 within a space A containing a set of objects described in an initial matrix 
M, comprising the construction of a subset co of objects contained in the said 
domain 8 by extracting a matrix jj. from the initial matrix M. 

According to the invention, the method comprises the following steps: 

- creation of a matrixing M of the space A by superimposing a number of 
geometric matrices with different specifications p (mesh sizes), representing 
coverage of the domain by a defined, homogeneous set of similar subdomains, 
each of the meshes of each geometric matrix being identified by a unique specific 
numerical value called the matrix code, 

- for the whole of the matrixing M, determination of all the meshes included in 
domain 8 or intersected by domain 8, and the number of relevant objects as the 
sum of the numbers of objects of the relevant meshes, 

- sorting of matrix M by matrix codes following a predetermined order, for 
example decreasing, of the specifications p, and 

- construction of the extraction matrix describing just the objects affected by the 
said meshes included in domain 8 or intersected by domain 8. 
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It should be noted that document US50301 17 describes a digital system 
for generating geographical maps, in which the cartographic data is organized 
according to a hierarchy of successive magnitudes or levels, in a pyramid shape 
with a small number of tiles (for example 4) at the top and 4 16 tiles at the base of 
5 a pyramid consisting for example of 16 levels. Even though diagrams in this 

document, in particular Figs. 8 and 13, show superimposed grids or matrices, it is 
not a question of a method of object retrieval according to a "multiple sieve" 
technique. 

In addition, from document US5694534, a device is known for storing a 
1 0 representation of topological structures and methods for constructing and 

searching this representation. The device described comprises a data storage 
medium, a database stored on this medium and containing blocks carrying data 
representing topological characteristics of the structure represented at a given 
level of detail. The data in each of these blocks are a representation of a carrier 
15 which is a closed set that includes a given topological object. In the method 
described, it is a matter of constructing connected chains of cells constituting 
topological subcomplexes. The concept of multiple sieving for purposes of 
constructing an extraction matrix is not disclosed in this document. 

In the search technique according to the invention, the objects can for 
20 example be sorted by increasing matrix codes. 

In an advantageous embodiment of the invention, the search technique 
additionally comprises storage of a list of matrix codes with a pointer to the first 
object of each code. 

Selection of the objects can be effected for example by a cursor which 
25 passes through the list of matrix codes. 

According to another aspect of the invention, a method of vectorial 
cartography is proposed, for mapping a territory A comprising a set of objects 
described by a matrix M; employing the method of object retrieval according to 
the invention, characterized in that it comprises the following steps: 



- determination of the list of active meshes, comprising calculation of the set of 
matrix codes corresponding to the meshes that intersect the search domain 8, 
and 

- selection of the objects and processing thereof if required. 

When this method is implemented in a two-dimensional domain, it then 
comprises matrixing according to a number of rectangular geometric matrices 
having regular meshes, each geometric matrix having a different specification 
including the length and the width of the meshes of the said matrix. 

For use in two-dimensional cartography extended to multiple levels, the 
method of vectorial cartography according to the invention then comprises the 
following sequences: 

- the initial set Q of N objects is broken down into a large number Z of 
subsets Clc, of same-level objects, with £ = 1 to Z, 

- the initial matrix M is broken down into Z matrices M ? , 

- matrixing is employed independently for each subset. 
According to yet another aspect of the invention, a system of vectorial 

cartography is proposed comprising means of processing objects contained in a 
domain 8 within a space A containing a set of objects described in an initial matrix 
M, comprising construction of a subset ra of objects contained in the said domain 
8 by extracting a matrix jj. from the initial matrix M, characterized in that it 
comprises in addition: 

- means for creating a matrixing M of the space A by superimposing a large 
number of geometric matrices with different specifications p, representing 
coverage of the domain with a defined and homogeneous set of similar 
subdomains, each of the meshes of each geometric matrix being identified by a 
unique specific index called the matrix code, 



- means of determining, for the whole of the matrixing M, all of the meshes 
included in domain 5 or intersected by domain 8, and the number of relevant 
objects as the sum of the number of objects with the relevant meshes, 

- means of sorting the matrix M by matrix codes according to a predetermined 
order, for example decreasing, of the specifications p, and 

- means of constructing the extraction matrix jj. describing just the objects 
affected by the said meshes included in the domain 8 or intersected by domain 8. 

Other advantages and characteristics of the invention will become clear on 
examining the detailed description of one way of implementing it, which is in no 
way limiting, and the appended drawings in which: 

Fig. 1 is a schematic representation of successive sieving carried out in the 

search method according to the invention; and 

Fig. 2 is a schematic diagram illustrating sorting of the matrix of objects from 
the set of objects by matrix codes 

The method proposed within the scope of the method of object retrieval 
according to the invention aims to speed up the process of obtaining the extract |_i 
without altering the objects in any way, so as not to increase their number and 
consequently the size of the matrix M. To achieve this, it takes account of the 
elementary space X occupied by each object, information that is contained in the 
matrix M, and uses it to organize successive "sieving", through "sieves" or 
"meshes" that become finer and finer. 

As in granulometry, groups of objects of particular sizes are obtained. 
Thus, referring to Fig. 1, if we consider a domain A containing objects 1, 2, 3, ...i, 
j...N, to which a set of grids of specifications p1, p2, p3,... pll is applied, the grid 
of specification p1 will "hold back" object 2, the next grid of specification p2 will 
hold back objects j and N, whereas the grid of specification p3 does not hold back 
any of the objects under consideration, with the last grid of specification pn finally 
collecting the object i and the point object 3. 



If the space X occupied by each object does not already exist in matrix M, 
it must be calculated. This space can, moreover, be recorded in matrix M, and in 
this case the size of M will increase by this information, but not the number of 
objects. 

For a given size of object, the "sieve" is a geometric matrix representing 
the coverage of domain A by a defined, homogeneous set of similar subdomains, 
not necessarily equal but close to a standardized occupied space p. Coverage of 
the domain does not have to be strict, in the sense that although the domain is 
effectively completely covered, there may nevertheless be intersections between 
the subdomains (between "meshes"). The geometric matrix representing the 
subset described is called a standard grid p. Each of the p of the meshes of the 
grid is identified by a unique specific numerical value called the matrix code of 
specification p, n p . By superimposing n grids with different specifications p, we 
create the matrixing M of depth n of domain A. This will comprise a total of 
meshes: 

^ = P 1 + -p2 + - + - P n (3.1) 

each one identified by its matrix code n pi j, with i = I, n (the depth, from I to n) 
and j = I, pj (the number of the mesh in the grid of depth i, from I to pi ). The first 
grid generally consists of a single mesh covering the whole domain and gathers 
together all the objects which, regardless of the space they occupy X, never "fall" 
into a finer mesh, but always "straddle" the whole sieve; then, of course, p i = I. 

Within a given matrixing M it is possible - and relatively easy - to 
determine, for each object of the set Q, the 



smallest mesh of specification p that contains it completely (the specification will 
always be greater than or equal to the space occupied: p > X) and hence ascribe 
to it a matrix code n pi j. Next we sort the matrix M of N objects of the set Q by 
matrix codes, for example in decreasing order of specifications p. We thus obtain 
an ordered description of the N objects in the sorted matrix M, ranging from those 
occupying the most space, to the "smallest" ones. For a given mesh, identified by 
its matrix code, there will be n objects; obviously, S n = N and certain n may be 
zero. The sort can be illustrated by the schematic diagram in Fig. 2. 

For construction of the subset co of o objects contained in domain 5, and 
therefore the extract \i from matrix M, we proceed in two steps: 

A. Determine, for the whole of the matrixing M (of depth n), all the meshes 
of all the grids included in and intersected by 5. For each grid (of 
specification pi) there will be pi meshes, pi < . pi from (3.1); the total number 
of grids involved will be 

¥ = P 1 + P 2 + ... + P n (3.2) 
and therefore \j/ < V F. The number of relevant objects, N, is the sum of the 
numbers of objects of the relevant meshes and therefore also N < N. The 
inequalities are greater the smaller 8 is, relative to A. 

B. With the matrix ordered as in the above scheme (or any other indexing 
system permitting direct access to the objects of a mesh, and therefore of 
a given matrix code), only the objects affected by the \\r meshes adopted 
are verified directly. The problem posed initially, i.e. determination of the 
subset o of the v objects contained in domain 8 and construction of the 
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matrix \i describing them, requires no more than examination of the N 
objects (and not the initial N). 

As in the case of tiling, the sort can be applied in each mesh. The first 
major advantage of matrixing relative to tiling is that none of the objects of Q is 
5 altered, and the size of the description matrix M therefore remains unchanged. 
The efficiency of matrixing does of course depend on the fineness of mesh, and 
therefore on the depth IT, just as the efficiency of tiling depends on I, but there is 
no other obstacle to increasing the depth. The only obvious condition to be 
observed is 

10 P n S p mi „>U (3.3) 

The choice of optimum depth is a function of the distribution of the spaces 
X occupied by the objects of Q and of the space A occupied by the latter, and 
hence of the field of application of the problem posed: 

in astrophysics A is much larger than A, ma x and therefore we can have a very 
15 fine mesh starting from the second grid; n will be small, even limited to 2; 

in cartography, X max = A and X, max » X m \n, while observing (3.3), n can be 
large and will be chosen in relation to the distribution of the Vs, all other 
constraints being the same. 

In the particular case of point objects, the objects occupying zero space {X 
20 = 0) are objects reduced to a point relative to domain A. The number of them 
varies depending on the field of application: very numerous in astrophysics, less 
numerous in cartography. Obviously, all of these point objects will be located in 
the meshes of the finest matrix, the grid of depth n. It should be noted that if Q 
contained such point objects exclusively, tiling and matrixing would be of 
25 equivalent efficiency. It will be noted that simple sorting combined with a 
dichotomizing search would also be very efficient. 
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We shall now give an example of implementation of the method of two- 
dimensional vectorial cartography. The territories (the domains A), for a given set 
Q of N objects, described by a matrix M, cover on average the extent of a region, 
or even a small country. The space A occupied by the domain is the diagonal of 
5 the rectangle [(X min , Y min ), (X ma x, Y max )] circumscribing it. The sides of the 
rectangle are respectively AX = X max - X min and AY = Y max - Y min . 

Each object is circumscribed by an "elementary" rectangle [(x m j n , y m in), 
(Xmax, ymax)], whose diagonal constitutes the space X occupied by the object. 
The matrixing M is made up of rectangular grids with regular meshes 
1 0 ensuring strict coverage. The specification pj of the grid of depth i is given by the 
length (on X) and the width (on Y) of the meshes: 

P ix = AX/> ix and p iy = AY/> iy (4.1) 
For the grid of depth i, we shall therefore have pi = > ix * >i y rectangular 
meshes. A series of n values for > ix and for > iy gives the matrixing M of depth n of 
15 domain A. 

In the example of implementation presented here, it has been stipulated, 
for reasons of conciseness, that all the matrix codes n pi j should be coded on a 
single octet. For simplicity, we chose an equal number of meshes in both 
directions > = >ix = >i y with the series of values 

20 ) = {1,2, 3,7, 13} (4.2) 

of depth n = 5 and which has, for the grid of the first depth, a single mesh equal 
to the territory. The use of a series of prime numbers guarantees that the 
boundaries of two meshes with different depths are never superimposed. Thus, 
an object whose occupied space is less than that of a mesh of depth i but whose 

25 geometry does not 
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allow it to be ascribed to a mesh at this depth (because it "straddles" at least two 
meshes of this depth), has every chance of being ascribed to a mesh at the next 
lower depth (i-1). Whatever the extraction domain, the number of inactive meshes 
(i.e. not intersecting the domain) increases with the depth (this increase is not 
necessarily strict, of course). The result is that the efficiency of a meshing system 
increases with increasing proportion of objects situated in the meshes with the 
greatest depths. 

The total number of meshes *F is 

1 2 + 2 2 + 3 2 + 7 2 + 13 2 = 232 (4.3) 
with the matrix codes indexable on an octet, from 0 to 231 . The mesh/index 
correspondence is as follows: 



depth 5: grid 13x13: 0 to 168 

depth 4: grid 7x7: 169 to 217 

depth 3: grid 3x3: 218 to 226 

depth 2: grid 2x2: 227 to 230 

depth 1: grid 1x1: 231 



Within the same grid, the meshes are sorted by increasing x and then by 
increasing y, as shown in Table 1 below. Extraction proper is made up of two 
steps: 

- determination of the list of active meshes, 

- selection of the objects and processing thereof if required. 

Determination of the active meshes merely consists of calculating the set 
of matrix codes corresponding to the meshes that intersect the search domain 5 
(in this case a rectangle). 

In the file, the objects are sorted by increasing matrix codes. A list of the 
matrix codes with a pointer to the first object of each code is also stored in the 
file. Selection of the objects is effected by a cursor that passes through the list of 
matrix codes. 
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Depth 5: 13x13 
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Depth 4: 7x7 
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When the cursor encounters an active code, the latter then passes through 
the set of objects with this code. Only the objects "seen" by the cursor are 
"processed". The processing itself can in its turn effect a sort, since the objects 
selected by the cursor do not necessarily intersect domain 8 (the fact that the 
mesh to which an object is associated intersects domain 5 is a condition that is 
necessary but not sufficient). This second sort is, naturally, heavier, since it 
requires reading of each object, then calculation of intersection with domain 8. 
We can thus understand the benefit of the preliminary sort effected by the cursor 
on the basis of the matrixing system. 

We shall now describe an example of application of the method of 
vectorial cartography with a pseudo-three-dimensional or two-dimensional extent 
(called "at levels"). 

In certain fields of application, such as "road-map" cartography (intended 
for applications in road-traffic guidance and navigation), it is not necessary to 
operate on a set of objects in 3D (with the Z for all the points), whose matrix M is 
inevitably more bulky than that of the set of the same objects represented on the 
flat (in 2D). In fact, to solve for example the problem of correct 2D graphical 
representation, it would be sufficient to know the relative vertical position of the 
objects: in short, it is sufficient to know who passes over whom (the "over/under" 
problem). 

The concept of level is then introduced: a level is made up of the set of 
objects located at one and the same relative vertical position. 

The method of vectorial cartography, when employed in an extended 
pseudo-three-dimensional or two-dimensional application, then comprises the 
following sequences: 

- The initial set Q of N objects is broken down into Z (total number of levels 
identified) subsets QL, of objects with the same level, with £ = 1 to Z 

- The initial matrix M is broken down into Z matrices M c . 
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- The matrixing described above is used independently for each subset. 
In other words, the initial set is first sorted by levels and then by spaces 

occupied, with the corresponding grids; in addition, this permits the use of 
different matrixing for different levels. 

Matrixing is employed in two separate stages in the method of searching 
for objects belonging to a subset: 

- a first stage comprising: 

- definition of the matrixing M itself, i.e. selection of the grids, 

- sorting of the objects of the set by increasing or decreasing order of their 
matrix code n p y, with i = 1 , n (the depth, from 1 to n) and j = 1 , . pi (the 
number of the mesh in the grid of depth i, from 1 to pi ). 

- a second stage comprising: 

- a step A of searching the "active" meshes (in this instance, the active 
matrix codes), i.e. those affected by 8, 

- for each active matrix code o p ij, a step B of searching the corresponding 
objects and of verification, for them alone, of membership of 8. 

We now present an example of implementation of matrix codes in the C++ 
language for the second stage of the method. So that the extracts of codes that 
follow can be understood, the definitions of four types and structures used are 
given below: 



// Types 

typedefTGF_INT32 
typedefTGF_UINT32 



long; 

unsigned long; 



// Structures 



typedef struct { 

TGFJNT32 



x; 



TGFJNT32 



y; 



}TGF_LONGPOINT; 
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typedef struct { 

TGF_LONGPOINT min; 

TGF_LONGPOINT max; 
} TGFLONGRECT; 

5 

The class CGFMatrixHandler contains all the methods that can be used in 
calculations connected with matrix codes. The method that follows permits 
determination of the list of "active" matrix codes (step A) relative to a domain (in 
this case a rectangle). The Boolean table mMatrixCodesList, class member data, 
10 is dimensioned by the total number of matrix codes. 

// Given inRect, compute list of "active" matrix codes 

bool * CGFMatrixHandler::GetMatrixCodesList(TGF_LONGRECT inRect) 

1 15 // Get global territory (domain) 

TGFJ.ONGRECT territory = Get TerritoryO; 

// Reset mMatrixCodesList. GetNbMatrixCodes returns total number of matrix 
%i codes memset(mMatrixCodesList, false, GetNbMatrixCodes()*sizeof(bool)); 

7 20 

M" // Now set "active " codes to "true" 

\ ^ TGF_LONGRECT currentMatrixSquare; 

TGF_UINT32 w, h, col, line, res, currentCode = 0; 

0 

hd-. 25 II GetCountQ returns matrix's "deepness" (number of levels) 

for (TGF_UINT32 i=0; i<GetCount(); i++) { 

// GetLevel(i) returns grid's resolution for level i 
res = GetLevel(i); 

30 // Compute matrix squares' dimension at level i 

w = (territory.max.x - territory.min.x + res-1) / res; 
h = ((territory. max.y - territory.min.y + res-1) / res; 

// initialize currentMatrixSquare 
35 currentMatrixSquare.min.x = territory, min. x; 

currentMatrixSquare. max.x = currentMatrixSquare.min.x+w; 



40 



// Loop on each matrix square of level i 
for (col=0; col<resol; col++){ 

currentMatrixSquare. min.y = territory.min.y; 
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currentMatrixSquare.max.y = currentMatrixSquare.min.y + h; 
for (line=0; line<resol; line++) { 

// Check if currentMatrixSquare intersects inRect. 

II If so, currentCode is "active" 

if (RectintersectsRect(currentMatrixSquare, inRect)) 
mMatrixCodesList[currentCode] = true; 



currentCode++; 

currentMatrixSquare.min.y += h; 
10 currentMatrixSquare.max.y += h: 

} 

currentMatrixSquare. min.x+= w; 
currentMatrixSquare. max. x+= w; 

} 

15 } 

return mMatrixCodesList; 

} 

The class CGFCursor permits access to the relevant objects (first part of 
20 step B), applying a filter to the matrix codes. The method Next() permits access 
to the next object in the file taking this filter into account. 



CGFCursor::CGFCursor(bool * inMatrixValues) 

: mMatrixCodes(inMatrixValues) 

25 { 

mCurMatrixCode = GetMatrixHandler()->GetNbMatrixCodes()-1; 
mCurMatrixCodelndex = (TGFJJINT32) -1; 
mCurMatrixCodeObjectsCount = (TGFJJINT32) -1; 
mCurObjectinMatrixCode = (TGF_UINT32) -1; 
30 mCurObjectOffset = (TGFJJINT32) -1 ; 

mCurObjectSize = 0; 
minitialized = false; 

} 

35 // Return offset of next object in file or (TGF_UINT32) -1 if none 

TGFJJINT32 CGFCursor::Next() 
{ 

for(;;){ 

// Move on to the next object 
40 for (;;) { 
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// Calculate the offset of this object. If current object's 
II offset is -1, then we are at the beginning of a matrix code 
II and we don't need to increment the previous offset to get to 
II the current one. 

5 if (mCurObjectOffset != (TGF_UINT32) -1) { 

// Next object in current matrix code 
if (++mCurObjectlnMatrixCode >= 

mCurMatrixCodeObjectsCount) break; 
// Compute the next object's position 
1 0 mCurObjectOffset += mCurObjectSize; 

} 

else{ 

mCurObjectlnMatrixCode = 0; 
mCurObjectOffset = 
1 5 GetFirstObjectlnNthMatrixCode(mCurMatrixCodelndex); 

} 

return mCurObjectOffset; 

} 

// Move on to the next matrix code 
20 NextMatrixCode(); 
} 

} 

bool CGFCursor::NextMatrixCode() 
{ 

25 // Reset the current offset so that the NextQ method 

II knows that it should start at the beginning of the matrix code. 
mCurObjectOffset = (TGF_UINT32) -1; 
mCurObjectSize = 0; 

mCurObjectlnMatrixCode = (TGFJJINT32) -1; 
30 mCurMatrixCodeObjectsCount = (TGFJJINT32) -1 ; 

// Handle the case where we are filtering on a matrix codes' list 
if (mMatrixCodes) { 
for(;;){ 

// Get the next code 

35 if (++mCurMatrixCodelndex >= GetMatrixCodesCount()) { 

mCurMatrixCodelndex = (TGF_UINT32) -1; 
mCurMatrixCode = GetMatrixHandler()-> 

GetNbMatrixCodesO -1; 
return false; 

40 } 

TGFJJINT32 code = GetNthMatrixCode(mCurMatrixCodelndex); 
// If this code is activated in the array, update variables 
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if (mMatrixCodesfcode}) { 

mCurMatrixCode = code; 
mCurMatrixCodeObjectsCount = 

5 GetObjectCountlnNthMatrixCode(mCurMatrixCodelndex); 

return true; 

} 

} 

} 

10 // Not filtering on matrix codes: get the next code 

if (++mCurMatrixCodelndex >= GetMatrixCodesCountO) { 
// Does not exist 

mCurMatrixCodelndex = (TGFJJINT32) -1 ; 
mCurMatrixCode = GetMatrixHandlerO ->GetNbMatrixCodes()-1 ; 
1 5 return false; 

} 

mCurMatrixCode = GetNthMatrixCode(mCurMatrixCodelndex); 
mCurMatrixCodeObjectsCount = 

GetObjectCountlnNthMatrixCode(mCurMatrixCodelndex); 

20 return true; 

} 

Access to the objects intersecting a rectangle rect (second part of step B) 
is then effected in the following way: 
{ 

25 TGF_LONGRECT rect, domain; 

TGF_UINT32 * matrixLevelList; 

// Set rect and get file's domain and matrixLevelList 

CGFMatrixHandler * mtxHdl = new(CGFMatrixHandler(domain, matrixLevelList)); 
30 CGFCursor * iter = new(CGFCursor(mtxHdl->GetMatrixCodesList(rect))); 

while (iter->NextO != (TGFJJINT32) -1) { 

// Read current object at offset iter->GetCurObjectOffsetQ 
II Check if object intersects rect (its matrix code does but not 
35 // necessarily the object itself)- 

II If so, process object 

} 

} 

A practical application of the search method according to the invention can 
40 be to carry out the first stage once and for all for 
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a given domain, before loading it into a reading device or equipment, especially if 
regular usage of the domain is anticipated. 

The second stage is always to be executed each time a search is 
requested and therefore permanent installation in the user terminal may be 
5 preferred. 

The extracts |i thus obtained by employing the method of vectorial 
cartography according to the invention can be stored in any type of data medium 
and incorporated in electronic equipment, such as equipment for communication 
and/or navigation. 

10 Of course, the invention is not limited to the examples that have just been 

described, and numerous adjustments can be made to these examples while 
remaining within the scope of the invention. 



